Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus.

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Comparison of Evolutionary Relationships between Branchiostoma floridae, Ciona intestinalis, and Homo sapiens Globins Provide Evidence of Gene Co-Option and Convergent Evolution.

Globins have been studied as model proteins to elucidate the principles of protein evolution. This was achieved by understanding the relationship between amino acid sequence, three-dimensional structure, physicochemical properties, and physiological function. Previous molecular phylogenies of chordate globin genes revealed the monophyletic evolution of urochordate globins and suggested convergent evolution. However, to provide evidence of convergent evolution, it is necessary to determine the physicochemical and functional similarities between vertebrates and urochordate globins. In this study, we determined the expression patterns of Ciona globin genes using real-time RT-PCR. Two genes (Gb-1 and Gb-2) were predominantly expressed in the branchial sac, heart, and hemocytes and were induced under hypoxia. Combined with the sequence analysis, our findings suggest that Gb-1/-2 correspond to vertebrate hemoglobin-α/-ß. However, we did not find a robust similarity between Gb-3, Gb-4, and vertebrate globins. These results suggested that, even though Ciona globins obtained their unique functions differently from vertebrate globins, the two of them shared some physicochemical features and physiological functions. Our findings offer a good example for understanding the molecular mechanisms underlying gene co-option and convergence, which could lead to evolutionary innovations.

Amphioxus as a model to study the evolution of development in chordates.

Cephalochordates and tunicates represent the only two groups of invertebrate chordates, and extant cephalochordates - commonly known as amphioxus or lancelets - are considered the best proxy for the chordate ancestor, from which they split around 520 million years ago. Amphioxus has been an important organism in the fields of zoology and embryology since the 18th century, and the morphological and genomic simplicity of cephalochordates (compared to vertebrates) makes amphioxus an attractive model for studying chordate biology at the cellular and molecular levels. Here we describe the life cycle of amphioxus, and discuss the natural histories and habitats of the different species of amphioxus. We also describe their use as laboratory animal models, and discuss the techniques that have been developed to study different aspects of amphioxus.

Amphioxus adenosine-to-inosine tRNA-editing enzyme that can perform C-to-U and A-to-I deamination of DNA.

Adenosine-to-inosine tRNA-editing enzyme has been identified for more than two decades, but the study on its DNA editing activity is rather scarce. We show that amphioxus (Branchiostoma japonicum) ADAT2 (BjADAT2) contains the active site 'HxE-PCxxC' and the key residues for target-base-binding, and amphioxus ADAT3 (BjADAT3) harbors both the N-terminal positively charged region and the C-terminal pseudo-catalytic domain important for recognition of substrates. The sequencing of BjADAT2-transformed Escherichia coli genome suggests that BjADAT2 has the potential to target E. coli DNA and can deaminate at TCG and GAA sites in the E. coli genome. Biochemical analyses further demonstrate that BjADAT2, in complex with BjADAT3, can perform A-to-I editing of tRNA and convert C-to-U and A-to-I deamination of DNA. We also show that BjADAT2 preferentially deaminates adenosines and cytidines in the loop of DNA hairpin structures of substrates, and BjADAT3 also affects the type of DNA substrate targeted by BjADAT2. Finally, we find that C89, N113, C148 and Y156 play critical roles in the DNA editing activity of BjADAT2. Collectively, our study indicates that BjADAT2/3 is the sole naturally occurring deaminase with both tRNA and DNA editing capacity identified so far in Metazoa.

Characterization of GSDME in amphioxus provides insights into the functional evolution of GSDM-mediated pyroptosis.

Members of the gasdermin (GSDM) family are pore-forming effectors that cause membrane permeabilization and pyroptosis, a lytic proinflammatory type of cell death. To reveal the functional evolution of GSDM-mediated pyroptosis at the transition from invertebrates to vertebrates, we conducted functional characterization of amphioxus GSDME (BbGSDME) and found that it can be cleaved by distinct caspase homologs, yielding the N253 and N304 termini with distinct functions. The N253 fragment binds to cell membrane, triggers pyroptosis, and inhibits bacterial growth, while the N304 performs negative regulation of N253-mediated cell death. Moreover, BbGSDME is associated with bacteria-induced tissue necrosis and transcriptionally regulated by BbIRF1/8 in amphioxus. Interestingly, several amino acids that are evolutionarily conserved were found to be important for the function of both BbGSDME and HsGSDME, shedding new lights on the functional regulation of GSDM-mediated inflammation.

Identification of a novel amphioxus leucine-rich repeat receptor involved in phagocytosis reveals a role for Slit2-N-type LRR in bacterial elimination.

The basal chordate amphioxus is a model for tracing the origin and evolution of vertebrate immunity. To explore the evolution of immunoreceptor signaling pathways, we searched the associated receptors of the amphioxus Branchiostoma belcheri (Bb) homolog of immunoreceptor signaling adaptor protein Grb2. Mass-spectrum analysis of BbGrb2 immunoprecipitates from B. belcheri intestine lysates revealed a folate receptor (FR) domain- and leucine-rich repeat (LRR)-containing protein (FrLRR). Sequence and structural analysis showed that FrLRR is a membrane protein with a predicted curved solenoid structure. The N-terminal Fr domain contains very few folate-binding sites; the following LRR region is a Slit2-type LRR, and a GPI-anchored site was predicted at the C-terminus. RT-PCR analysis showed FrLRR is a transcription-mediated fusion gene of BbFR-like and BbSlit2-N-like genes. Genomic DNA structure analysis implied the B. belcheri FrLRR gene locus and the corresponding locus in Branchiostoma floridae might be generated by exon shuffling of a Slit2-N-like gene into an FR gene. RT-qPCR, immunostaining, and immunoblot results showed that FrLRR was primarily distributed in B. belcheri intestinal tissue. We further demonstrated that FrLRR localized to the cell membrane and lysosomes. Functionally, FrLRR mediated and promoted bacteria-binding and phagocytosis, and FrLRR antibody blocking or Grb2 knockdown inhibited FrLRR-mediated phagocytosis. Interestingly, we found that human Slit2-N (hSlit2-N) also mediated direct bacteria-binding and phagocytosis which was inhibited by Slit2-N antibody blocking or Grb2 knockdown. Together, these results indicate FrLRR and hSlit2-N may function as phagocytotic-receptors to promote phagocytosis through Grb2, implying the Slit2-N-type-LRR-containing proteins play a role in bacterial binding and elimination.

SUMO and PIAS repress NF-κB activation in a basal chordate.

Small ubiquitin-like modifier (SUMO) regulates various biological processes, including the MyD88/TICAMs-IRAKs-TRAF6-NF-κB pathway, one of the core immune pathways. However, its functions are inconsistent between invertebrates and vertebrates and have rarely been investigated in lower chordates, including amphioxus and fishes. Here, we investigated the SUMOylation gene system in the amphioxus, a living basal chordate. We found that amphioxus has a SUMOylation system that has a complete set of genes and preserves several ancestral traits. We proceeded to study their molecular functions using the mammal cell lines. Both amphioxus SUMO1 and SUMO2 were shown to be able to attach to NF-κB Rel and to inhibit NF-κB activation by 50-75% in a dose-dependent fashion. The inhibition by SUMO2 could be further enhanced by the addition of the SUMO E2 ligase UBC9. In comparison, while human SUMO2 inhibited RelA, human SUMO1 slightly activated RelA. We also showed that, similar to human PIAS1-4, amphioxus PIAS could serve as a SUMO E3 ligase and promote its self-SUMOylation. This suggests that amphioxus PIAS is functionally compatible in human cells. Moreover, we showed that amphioxus PIAS is not only able to inhibit NF-κB activation induced by MyD88, TICAM-like, TRAF6 and IRAK4 but also able to suppress NF-κB Rel completely in the presence of SUMO1/2 in a dose-insensitive manner. This suggests that PIAS could effectively block Rel by promoting Rel SUMOylation. In comparison, in humans, only PIAS3, but not PIAS1/2/4, has been reported to promote NF-κB SUMOylation. Taken together, the findings from amphioxus, together with those from mammals and other species, not only offer insights into the functional volatility of the animal SUMO system, but also shed light on its evolutionary transitions from amphioxus to fish, and ultimately to humans.

Genomic Analysis of Amphioxus Reveals a Wide Range of Fragments Homologous to Viral Sequences.

Amphioxus species are considered living fossils and are important in the evolutionary study of chordates and vertebrates. To explore viral homologous sequences, a high-quality annotated genome of the Beihai amphioxus (Branchiostoma belcheri beihai) was examined using virus sequence queries. In this study, 347 homologous fragments (HFs) of viruses were identified in the genome of B. belcheri beihai, of which most were observed on 21 genome assembly scaffolds. HFs were preferentially located within protein-coding genes, particularly in their CDS regions and promoters. A range of amphioxus genes with a high frequency of HFs is proposed, including histone-related genes that are homologous to the Histone or Histone H2B domains of viruses. Together, this comprehensive analysis of viral HFs provides insights into the neglected role of viral integration in the evolution of amphioxus.

Retinoic Acid and POU Genes in Developing Amphioxus: A Focus on Neural Development.

POU genes are a family of evolutionarily conserved transcription factors with key functions in cell type specification and neurogenesis. In vitro experiments have indicated that the expression of some POU genes is controlled by the intercellular signaling molecule retinoic acid (RA). In this work, we aimed to characterize the roles of RA signaling in the regulation of POU genes in vivo. To do so, we studied POU genes during the development of the cephalochordate amphioxus, an animal model crucial for understanding the evolutionary origins of vertebrates. The expression patterns of amphioxus POU genes were assessed at different developmental stages by chromogenic in situ hybridization and hybridization chain reaction. Expression was further assessed in embryos subjected to pharmacological manipulation of endogenous RA signaling activity. In addition to a detailed description of the effects of these treatments on amphioxus POU gene expression, our survey included the first description of Pou2 and Pou6 expression in amphioxus embryos. We found that Pit-1, Pou2, Pou3l, and Pou6 expression are not affected by alterations of endogenous RA signaling levels. In contrast, our experiments indicated that Brn1/2/4 and Pou4 expression are regulated by RA signaling in the endoderm and the nerve cord, respectively. The effects of the treatments on Pou4 expression in the nerve cord revealed that, in developing amphioxus, RA signaling plays a dual role by (1) providing anteroposterior patterning information to neural cells and (2) specifying neural cell types. This finding is coherent with a terminal selector function of Pou4 for GABAergic neurons in amphioxus and represents the first description of RA-induced changes in POU gene expression in vivo.

Identification of Nodal-dependent enhancer of amphioxus Chordin sufficient to drive gene expression into the chordate dorsal organizer.

The core molecular mechanisms of dorsal organizer formation during gastrulation are highly conserved within the chordate lineage. One of the key characteristics is that Nodal signaling is required for the organizer-specific gene expression. This feature appears to be ancestral, as evidenced by the presence in the most basally divergent chordate amphioxus. To provide a better understanding of the evolution of organizer-specific gene regulation in chordates, we analyzed the cis-regulatory sequence of amphioxus Chordin in the context of the vertebrate embryo. First, we generated stable zebrafish transgenic lines, and by using light-sheet fluorescent microscopy, characterized in detail the expression pattern of GFP driven by the cis-regulatory sequences of amphioxus Chordin. Next, we performed a 5'deletion analysis and identified an enhancer sufficient to drive the expression of the reporter gene into a chordate dorsal organizer. Finally, we found that the identified enhancer element strongly depends on Nodal signaling, which is consistent with the well-established role of this pathway in the regulation of the expression of dorsal organizer-specific genes across chordates. The enhancer identified in our study may represent a suitable simple system to study the interplay of the evolutionarily conserved regulatory mechanisms operating during early chordate development.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications.

BACKGROUND: Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates. RESULTS: We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements. CONCLUSIONS: In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Localization of trypsin, chymotrypsin and elastase in the digestive tract of amphioxus Branchiostoma japonicum with implications to the origin of vertebrate pancreas.

The organ pancreas is characteristic of all vertebrates, but its evolutionary origin remains largely unknown. The serine proteinases trypsin, chymotrypsin and elastase are produced primarily in the pancreas in vertebrates. Here we clearly demonstrate by enzymatic activity assay, histochemical staining and in situ hybridization that trypsin, chymotrypsin and elastase are widely distributed in the digestive tract of the amphioxus Branchiostoma japonicum, especially in the hepatic caecum and mid-gut. This suggests that an extensive region of the amphioxus digestive tract including the hepatic caecum is homologous to the vertebrate exocrine pancreatic cells, providing a new angle for the study of the origin and evolution of vertebrate pancreas.

Highly distinct genetic programs for peripheral nervous system formation in chordates.

BACKGROUND: Vertebrates develop their peripheral nervous system (PNS) from transient unique embryonic structures, the neural crest, and the ectodermal placodes that are located at the border of the forming central nervous system. By contrast, in the invertebrate chordates, amphioxus and ascidians, a large part of the PNS originates at the opposite of the embryo, in the ventral ectoderm. In both groups, a biphasic mechanism regulates ventral PNS formation: high BMP levels specify a neurogenic territory within which glutamatergic epidermal sensory neuron formation is controlled by the Notch pathway. Given these similarities and the phylogenetic relationships within chordates, it is likely that ventral PNS is an ancestral feature in chordates and that it has been lost in vertebrates. RESULTS: In order to get insights into the molecular control of ventral PNS formation and to test the hypothesis of their homology and potential contribution to the emergence of vertebrate PNS, we undertook a close comparison of ventral PNS formation in the ascidian Phallusia mammillata and the amphioxus Branchiostoma lanceolatum. Using timed RNA-seq series, we identified novel markers of the ventral PNS during different phases of its development in both species. By extensively determining the expression of paralogous and orthologous genes, we observed that only a minority of genes have a shared expression in the ventral PNS. However, a large fraction of ventral PNS orthologous genes are expressed in the dorsally forming PNS of vertebrates. CONCLUSIONS: Our work has significantly increased the molecular characterization of ventral PNS formation in invertebrate chordates. The low observed conservation of gene expression in the ventral PNS suggests that the amphioxus and ascidian ventral PNS are either not homologous, or alternatively extensive drift has occurred in their regulatory mechanisms following a long period (600 My) of separate evolution and accelerated evolution in the ascidian lineage. The homology to genes expressed in the dorsally forming PNS of vertebrates suggests that ancestral sensory neurons gene networks have been redeployed in vertebrates.

Joint profiling of gene expression and chromatin accessibility during amphioxus development at single-cell resolution.

Vertebrate evolution was accompanied by two rounds of whole-genome duplication followed by functional divergence in terms of regulatory circuits and gene expression patterns. As a basal and slow-evolving chordate species, amphioxus is an ideal paradigm for exploring the origin and evolution of vertebrates. Single-cell sequencing has been widely used to construct the developmental cell atlas of several representative species of vertebrates (human, mouse, zebrafish, and frog) and tunicates (sea squirts). Here, we perform single-nucleus RNA sequencing (snRNA-seq) and single-cell assay for transposase accessible chromatin sequencing (scATAC-seq) for different stages of amphioxus (covering embryogenesis and adult tissues). With the datasets generated, we constructed a developmental tree for amphioxus cell fate commitment and lineage specification and characterize the underlying key regulators and genetic regulatory networks. The data are publicly available on the online platform AmphioxusAtlas.

PACAP/GCGa Is an Important Modulator of the Amphioxus CNS-Hatschek's Pit Axis, the Homolog of the Vertebrate Hypothalamic-Pituitary Axis in the Basal Chordates.

The Hatschek's pit in the cephalochordate amphioxus, an invertebrate deuterostome basal to chordates is suggested to be the functional homolog structure of the vertebrate adenohypophysis based on anatomy and expression of homologous neuroendocrine genes. However, the endocrine potential of the cephalochordate Hatschek's pit remains to be demonstrated as well as the physiological actions of the secreted neuropeptides. In this study, we have explored the distribution and characterize the potential function of the amphioxus PACAP/GCG precursor, which is the ortholog of the hypothalamic PACAP neuropeptide in vertebrates. In amphioxi, two PACAP/GCG transcripts PACAP/GCGa and PACAP/GCGbc that are alternative isoforms of a single gene with different peptide coding potentials were isolated. Immunofluorescence staining detected their expression around the nucleus of Rohde, supporting that this structure may be homologous of the neurosecretory cells of the vertebrate hypothalamus where abundant PACAP is found. PACAP/GCGa was also detected in the infundibulum-like downgrowth approaching the Hatschek's pit, indicating diffusion of PACAP/GCGa from the CNS to the pit via the infundibulum-like downgrowth. Under a high salinity challenge, PACAP/GCGa was upregulated in amphioxi head and PACAP/GCGa treatment increased expression of GHl in Hatschek's pit in a dose-dependent manner, suggesting that PACAP/GCGa may be involved in the regulation of GHl via hypothalamic-pituitary (HP)-like axis similar as in the vertebrates. Our results support that the amphioxus Hatschek's pit is likely to be the functional homolog of pituitary gland in vertebrates.

Functions of the FGF signalling pathway in cephalochordates provide insight into the evolution of the prechordal plate.

The fibroblast growth factor (FGF) signalling pathway plays various roles during vertebrate embryogenesis, from mesoderm formation to brain patterning. This diversity of functions relies on the fact that vertebrates possess the largest FGF gene complement among metazoans. In the cephalochordate amphioxus, which belongs to the chordate clade together with vertebrates and tunicates, we have previously shown that the main role of FGF during early development is the control of rostral somite formation. Inhibition of this signalling pathway induces the loss of these structures, resulting in an embryo without anterior segmented mesoderm, as in the vertebrate head. Here, by combining several approaches, we show that the anterior presumptive paraxial mesoderm cells acquire an anterior axial fate when FGF signal is inhibited and that they are later incorporated in the anterior notochord. Our analysis of notochord formation in wild type and in embryos in which FGF signalling is inhibited also reveals that amphioxus anterior notochord presents transient prechordal plate features. Altogether, our results give insight into how changes in FGF functions during chordate evolution might have participated to the emergence of the complex vertebrate head.

Deep whole-genome resequencing sheds light on the distribution and effect of amphioxus SNPs.

BACKGROUND: Amphioxus is a model organism for vertebrate evolutionary research. The significant contrast between morphological phenotypic similarity and high-level genetic polymorphism among amphioxus populations has aroused scientists' attention. Here we resequenced 21 amphioxus genomes to over 100X depth and mapped them to a haploid reference. RESULTS: More than 11.5 million common SNPs were detected in the amphioxus population, which mainly affect genes enriched in ion transport, signal transduction and cell adhesion, while protein structure analysis via AlphaFold2 revealed that these SNPs fail to bring effective structural variants. CONCLUSIONS: Our work provides explanation for "amphioxus polymorphism paradox" in a micro view, and generates an enhanced genomic dataset for amphioxus research.

The Origination of Growth Hormone/Insulin-Like Growth Factor System: A Story From Ancient Basal Chordate Amphioxus.

The growth hormone/insulin-like growth factor (GH/IGF) system, also called the pituitary-liver axis, has a somatotrophic role in the body. Although the GH/IGF system has always been regarded as a vertebrate-specific endocrine system, its actual origin remained unknown for a long time. The basal chordate, amphioxus, occupies an evolutionary position between vertebrates and invertebrates. Impressively, most of the members of the GH/IGF system are present in the amphioxus. The GH-like molecule in the amphioxus is mainly expressed in Hatschek's pit. It functions similarly to vertebrate GH and has a GH receptor-like binding partner. The amphioxus IGF-like peptide shows mitogenic activity and an expression pattern resembling that of vertebrate IGF-I. The receptor of IGF-like peptide and IGF binding protein (IGFBP) have also been demonstrated to exist in the amphioxus. These results reveal the origin of the gene families in the GH/IGF system, providing strong evidence that this system emerged in the amphioxus.

Gain of gene regulatory network interconnectivity at the origin of vertebrates.

SignificanceIn this manuscript, we address an essential question in developmental and evolutionary biology: How have changes in gene regulatory networks contributed to the invertebrate-to-vertebrate transition? To address this issue, we perturbed four signaling pathways critical for body plan formation in the cephalochordate amphioxus and in zebrafish and compared the effects of such perturbations on gene expression and gene regulation in both species. Our data reveal that many developmental genes have gained response to these signaling pathways in the vertebrate lineage. Moreover, we show that the interconnectivity between these pathways is much higher in zebrafish than in amphioxus. We conclude that this increased signaling pathway complexity likely contributed to vertebrate morphological novelties during evolution.

The invertebrate chordate amphioxus gives clues to vertebrate origins.

Amphioxus (cepholochordates) have long been used to infer how the vertebrates evolved from their invertebrate ancestors. However, some of the body part homologies between amphioxus and vertebrates have been controversial. This is not surprising as the amphioxus and vertebrate lineages separated half a billion years ago-plenty of time for independent loss and independent gain of features. The development of new techniques in the late 20th and early 21st centuries including transmission electron microscopy and serial blockface scanning electron microscopy in combination with in situ hybridization and immunocytochemistry to reveal spatio-temporal patterns of gene expression and gene products have greatly strengthened inference of some homologies (like those between regions of the central nervous system), although others (like nephridia) still need further support. These major advances in establishing homologies between amphioxus and vertebrates, together with strong support from comparative genomics, have firmly established amphioxus as a stand-in or model for the ancestral vertebrate.